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Entangled states are the basic ingredients of quantum information and quantum com- 
munication. For pure states of any dimension and any number of particles it is easy to 
determine if a state is entangled or not, the property is revealed by its reduced density 
matrices (see, e.g., Ref. For mixed states, however, it is in general difficult to decide 
whether a given density matrix is entangled or separable. Since quantum mechanical states 
are in practice rather mixed than pure, one seeks for criteria to detect entanglement and 
separability. Additionally, entangled states can be classified according to their distillability 
- a state is called distillahle if from many copies of the same not maximally entangled state 
we can distill fewer maximally entangled states with statistical local operations and classical 
communication (SLOCC). Interestingly, not all entangled states are distillable. For states 
of higher dimension than 2x2 and 2x3, there exist states that are entangled but not dis- 
tillable, so-called bound entangled states 

BiHiHB- Apart from being a mathematical 
oddity bound entangled states physically characterize a set of states that are an obstacle for 
quantum information processing tasks, since once a system evolved into a bound entangled 
state, there is no way to retain a maximally entangled state with SLOCC. Nevertheless, 
concerning their usefulness, bound entangled states are not equivalent to separable states 
since their entanglement can be "activated" when used together with distillable entangled 
states 0]. It is the intention of this paper to present a certain family of two-qutrit states 
that lie within the magic simplex jl] - this family reveals a high geometric symmetry, can 
be classified according to its entanglement properties and contains large regions of bound 
entanglement. 

We consider a Hilbert-Schmidt space Aa ® Ab of operators on the finite dimensional 
bipartite Hilbert space HA^Ti-B, with dimension (Ia x (Ib, D := dAds- States p (i.e. density 
matrices) are elements of Aa®Ab with the properties = p, Tr p = 1 and p > 0. A scalar 
product on Aa ® Ab is defined by {A, B) = TrA^B with A,BeAA® Ab- 

Two basic concepts in the context of entanglement detection are entanglement witnesses 
[sl, 10, [m, [H, 13] and the PPT Criterion (positive partial transposition) jl, 14|. An en- 



tanglement witness A detects the entanglement of states due to the convexity of the set of 
separable states S: A state p is entangled iff there exists a Hermitian operator A such that 

{p,A) = TipA < 0, (1) 
{a,A) = TTaA > \/ a e S . (2) 
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For a given state p it is easy to find operators that provide Eq. ([T]), but in general difficult to 
show Eq. ([2]). Nevertheless it proves useful in many cases. If there exists a separable state a 
for which Tr a A = 0, then A is called optimal entanglement witness. Optimal entanglement 
witnesses correspond to hyperplanes given by TraA = which are tangent to the set of 
separable states S. 

An operational criterion that is easy to apply is the PPT Criterion: A separable state a 
stays positive under partial transposition, a^'^ := (1 T) a > . We call a state PPT that 
is positive under partial transposition, and NPT a state that is not. Note that the PPT 
criterion is a necessary criterion for separability [13], that means if a state p is NPT, it has 
to be entangled. But if it is PPT, this does not automatically imply that it is separable, this 
is true for dimensions 2x2 and 2x3 only j^. Amazingly, entangled PPT states are not 
distillable, thus bound entangled jl, H, 0] ■ Clearly the PPT Criterion cannot distinguish PPT 
entangled from separable states, therefore other entanglement criteria have to be used, e.g. 
the entanglement witness criterion, which is the method that we will apply in this paper. 

In Ref. we presented a method to detect bound entanglement that is based on a 
geometrical approach: Consider a line of states p\ inside the convex set of the PPT states 
between a PPT state pppr and the maximally mixed state. 

Pa := Apppt + ^^^^Id, 0<A<1, (3) 

and an operator 

Ca = Pa - Pppt - (pa, Pa - Pppt) Id • (4) 

The operator Cx is constructed geometrically in such way that we have Tr C\px = and 
Tr CaPppt < (here A < 1, the limit A — > 1 is discussed later on). It divides the whole 
state space into states p„ for which Tr Cxpn < and states pp for which Tr CxPp > 0, see 



Ref. [15|. If Cx is an entanglement witness then all states pn on one side of the plane are 
entangled. For a particular entanglement witness Cj^ all states px are PPT and entangled 
for A < A < 1, thus bound entangled. Clearly the difficulty of the presented method lies 
in proving that Ca > V a G 5*. Different hues (13]) are distinguished by different "starting 
states" Pppt which will be called "starting points" throughout this text. 

Next we apply the above method to a three-parameter family of the magic simplex of 
two qutrits. The magic simplex W of bipartite qutrits (dimension 3 x 3) is introduced in 
Ref. [l6|, extended to qudits (dimension dxd) in Ref. [8|, reviewed and discussed in Ref. [^. 
It is defined as the set of all states that are a mixture of Bell states Pr, 



nmi 



W := < ^ qnmPnm \ qnm > 0, ^ ^nm = 1 j" , (5) 
[^n,m=0 n,m ) 

where 

Pn^m := (f/„™®l)|0+)(0+|(f/L®l)- (6) 

The vector state denotes the maximally entangled state = J2j |j) ® Ij) 
the Weyl operators Unm = Ylk=o '^'^ \k) {(k + m) mod d\ form an orthogonal basis {Unm} 



{n,m = 0,1, d — 1) of the Hilbert-Schmidt space [17|, [18 

The indices n and m can be viewed as coordinates in a discrete phase space, and the points 
for fixed n or m constitute lines within. The magic simplex is hi ghly symmetric and exhibits 
the same geometry under a line change (for details see Refs. [^, ly]). It thus represents the 
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eight-dimensional analogue of the three-dimensional simplex, a tetrahedron, of two-qubits 
which is the set of states that are mixtures of the four two-qubit Bell states |12l . Il9l . 



The two-qudit Bell states ([6]) play an important role in extensions of quantum communica- 
tion procedures to higher dimensional systems - e.g. in the quantum teleportation protocol 
13] . They constitute the corresponding maximally entangled orthogonal basis and the Weyl 



operators f/„m are a generalization of the Pauli operators, in the sense that Bob has to apply 
one of the cP operators (including the identity operator) to obtain the teleported state. 

We want to consider a three-parameter family of states, introduced in Ref. [l3|, of the 
magic simplex that are a mixture of the maximally mixed state and two phase space lines, 

1 — <y — 3 — 'y 3 i 

PaAi ■■= ^ ^1 + + ^ (Pio + P20) + j (Poi + Pu + P21) ■ (7) 

The parameters are constrained by the positivity requirement Pa,/?,^ > 0, which results in 
a < 1(3/2 + 1-7, a < -/3 + 1 - 7, a < -/3 + 1 + 27, and a > /3/8 - 1/8 + 7/8. 
The constraints geometrically represent a pyramid with triangular base, see Fig. [H The 
PPT Criterion selects those points (a, /?, 7) that correspond to positive operators when 
subjected to partial transposition. We obtain the constraints a < —(3 — 1/2 + 7/2, as well 
as a < (-2 + 11/? - 7 + 3yA)/16 and a< (-2 + 11/3 - 7 - 3v^)/16, where we defined 
A = 4 + 9/5^ + 47 — 77^ — 6/3(2 + 7), which form a cone that intersects the pyramid, and in 
the intersection region lie the PPT states, see Fig. [H 

A parameterized line of states Pa,/3,7 situated on the boundary plane a = 7/3/2 + 1 — 7 
(which we will from now on refer to as the boundary plane) is given by 

6 - 6 ^ 2b 5-26 
pb = Pa,/3,7 with a = 13 = - — , 7 = — — , (8) 

which exactly corresponds to the states 

P6 = ^|0'+>(0'+| + + 0<6<5 (9) 

that were introduced in Ref. [3]. We call this one-parameter family Horodecki states or 
Horodecki line. It is shown in Ref. that these states are entangled for 4 < 6 < 5, bound 
entangled for 3 < 6 < 4 and separable for 2 < 6 < 3. For convenience we express the 
parameter b in terms of 7 to obtain a parametrization of the Horodecki line with 7. 



In Ref. [13| we discovered planar regions of bound entanglement that go beyond the 
Horodecki line with the method discussed before, where the starting points pppT of the line 
Pa (E]) are chosen as the bound entangled Horodecki states. The inequality ([2]) is proven 
by expressing the operator Cx and the separable states in terms of Weyl operator tensor 
products. In this way we arrive at a sufficient condition for Cx such that Eq. ([2]) is satisfied: 

Lemma 1. For any Hermitian operator C of a bipartite Hilbert- Schmidt space of dimension 
d X d that is of the form 

C = and -1)1^2 + ^ Cnm Unm ® U^nm j , « ^ R+, Cnm ^ C (lO) 
\ n,m=0 / 

the expectation value for all separable states is positive, 

{PiC) > Q WpES, if \cnm\ < 1 Wn,m. (11) 
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FIG. 1: (Color online) In the parameter space representing the Hilbert-Schmidt space geometry 
the states /Oa,/3,7 © he inside a pyramid with triangular base, the PPT points form a cone. The 
intersection region includes the PPT states, i.e. separable and bound entangled states. The 
maximally mixed state is represented by the black dot. The line (blue) on the boundary plane 
indicates the Horodecki states ([9]), the planar section (red) jutting out of the Horodecki line shows 
the bound entangled states that are detected using starting points on the Horodecki line only. The 
protruding volume (blue/green) is the volume of bound entangled states that can be obtained with 
starting points deviating from the Horodecki line. In fact the picture is symmetrical for positive 
and negative values of 7, but for the sake of clarity only one volume is depicted. All parameter 
axes are chosen non-orthogonal to each other such that they become orthogonal to the boundary 
of the positivity region in order to reproduce the symmetry of the magic simplex. 

If we choose the starting points not only on the Horodecki line, but also in a neighborhood 
of the line on the boundary plane, we can detect volumes of bound entanglement, see Fig. [H 
These starting points Ppiane are parameterized with an additional parameter e that accounts 
for the amount of deviation from the line. 



The parameters e and 7 of the states Ppiane are restricted by the condition \cnm\ < 1 ^n,m 
of Lemma [H The coefficients Cnmil) of C\ fulfil Cnmil) = c*m(— 7), hence Lemma [1] can be 
applied for positive and negative values of 7, and we obtain a symmetric picture. We restrict 
ourselves to positive values of 7, obtaining the "negative side" by applying 7 —7. The 
starting points on the plane can be seen in Fig. [T] as the points where the volumes of bound 
entangled states emerge from the boundary plane. Note that only the states inside the 
volumes are bound entangled for sure, states p^^^^^^ on the surface can be either separable 
or bound entangled. The reason is that for points p^"''f^'=*^ we have Tr p^urface^surface _ q 
and in general we do not know whether the operator C^'^^^'^'^ defines a tangent plane to the 



Pplane — Pa,f3,"/ 



with 




e e R. 



(12) 
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FIG. 2: (Color online) The entanglement properties of states Pa,i3,'y on the boundary plane a = 
(7/2)/3 + 1 — 7. The separable states are situated in the triangular region (red), the PPT entangled 
(bound entangled) states lie in between the triangular and parabolic region (orange), and the 
remaining states are NPT entangled. Interestingly, the Horodecki line (blue) runs through all 
entanglement characteristics. 

set of separable states, i.e. whether it is already optimal. The detected regions of bound 
entanglement, the "old" planar region and the "new" volumes, are plotted in Fig. [TJ 

Entanglement witnesses of the form (j4]) correspond to planes in our three-dimensional 
illustration of the Hilbert-Schmidt space. They detect not only the entanglement of states 
on lines p\ ([3]), but also of all states lying on the same side of the plane that includes 
the starting points pppT of px- Consider for example a starting point Ppi^ne on the 
boundary plane that represents the lower tip of the bound entangled volume in Fig. [1], 
it is given by (e = —1/4,7 = 1/4)- Using Lemma [1] we find that for Ppiane the operator 
Ci := Cx/X{l — A) is an entanglement witness for A ^ 1 and has a simple matrix form 
which we can easily calculate from Eq. (jl]). 



Now we explain how to detect further bound entangled states and the region of separabil- 
ity. The trace TiCiPp^^^^ = provides the plane equation Pli : a = 2(1 + 2/3 — 7)/5. Plane 
Pli intersects the boundary plane on the line la ■ P = 2(— 1 + 7)/9, the cone of PPT states 
intersects the boundary plane on the curve /ft : /3 = (-4 + 37 + v^4- 872)79, and and 
lb cross each other at 7 = and 7 = 1. These crossing points are separable states: 7 = 
corresponds to two-parameter states 13, 1^, in this case it is proved in Ref. [16] that all 



PPT states are separable, and for 7 = 1 we have a = (3 = 0, i.e. an equal mixture of the Bell 



states Poi, Pu, and P21, which is separable too [16|. The hyperplane of the entanglement 
witness cuts the boundary plane at the line la, it contains two separable states, therefore 
it represents an optimal entanglement witness that has to be tangent to the set of separa- 
ble states S. Since S is convex the line has to describe the boundary of the separable 
states between 7 = and 7 = 1. Therefore the whole region between la and has to be 
bound entangled, it is given by < 7 < 1, 2(-l + 7)/9 < /3 < (-4 + 87 + a/4 - 8779) 
and a = 7/?/2 + l — 7. For the Horodecki states we thus find bound entanglement for 
1/7 < 7 < 8/7 and separability for < 7 < 1/7, which (for 7 — —7) recovers the result 
of Ref. U The regions of separable, PPT entangled and NPT entangled states on the 
boundary plane are illustrated in Fig. [2l 

The plane Pli representing the entanglement witness also intersects the PPT cone inside 
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FIG. 3: (Color online) Right: The detected regions of bound entanglement extending from the 
boundary plane (enlarged) as intersections of the planes Pli (green), PI2 (red) and P/3 (blue) 
(from bottom to top) corresponding to three entanglement witnesses. Left: The polygon (green) 
encloses the states that are necessarily separable. In between the polygon and the surface of the 
PPT cone lie the detected volumes of bound entanglement. 

the pyramid; thus the bound entangled states extend from the boundary plane into the 
volume. The region of bound entanglement is then confined by the surface of the PPT cone 
a = (—2 + lip — 7 + 3v^)/16, the plane Ph, and the boundary plane. 

We can detect regions of bound entanglement that reach even further into the pyra- 
mid by using other entanglement witnesses that correspond to other planes intersect- 
ing the cone. We construct again witnesses with our geometric method on lines start- 
ing on the boundary plane around the Horodecki line and such we obtain the small- 
est possible value of A, i.e., = (3 + vT3)/8 ~ 0.826, for the starting point with 
eo = (-25 + 7v^)/2 ~ 0.12 and 7 = v^S + lleo/3 - 5eg/12 / 7 ~ 0.35. The state 
on this line with AJ:^*^ is the outermost point of the bound entangled volume (the green 
tip in Fig. [1]) that reaches into the PPT cone. The entanglement witness C\ correspond- 
ing to A*^*n gives the plane equation Pk : a = (-524 + 148v^)-^[16 (-7 + 2v^) - 
4 (-5 + \/T3) (3+ (-94 + 26VT3 + 3(-5 + ^13)^/2^) 7]. It offers a new boundary for the 
bound entangled region that extends further into the pyramid (but detects fewer bound 
entangled states on the boundary plane of the pyramid). A final extension we find by con- 
sidering the line that starts at a point 7 = 2/7 of the intersection curve of Pli with the 
PPT cone. The minimal A such that Cx (jl]) represents an entanglement witness accord- 
ing to Lemma [1] is Amm = 7(2328 + 331v^ )/32763 ~ 0.94, and the corresponding plane 
Ph : a = (150 - 18v^)-^ [24 - 2V39 + (-42 + 9^39) /5 - 6 (-5 + ^39) 7] provides a 
new boundary for the bound entangled states reaching yet a bit further into the volume of 
the pyramid. The planes Pli, PI2, PI3 together with the detected volumes of bound entan- 
gled states we have depicted in Fig. [31 

Naturally the question arises if we can detect all regions of NPT entangled, PPT entangled 
and separable states for our three-parameter family of states ((71). As a first step to an 
answer we can identify a polygon of states that have to be separable for sure by connecting 
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all outermost separable states we know so far in a convex manner. These known states are 
the boundary of the PPT states for 7 = 0, which form a trapezoid, and the boundary of the 
separable states on the boundary plane. The resulting polygon can be seen in Fig. [31 

Numerical calculations imply that the set of separable states is larger than the constructed 
polygon, and that bound entanglement is mainly concentrated in the region we detected. 
These results require numerical minimizations to obtain the needed entanglement witnesses, 
and the detailed procedure will be discussed in a forthcoming article. 

Summarizing, we find new, large regions of bound entangled states in a three-parameter 
family of the magic simplex of two qutrits. We detect the definite regions of separable, PPT 
entangled and NPT entangled states on a boundary plane of the family of states by using 
entanglement witnesses that are constructed geometrically. The regions of bound entangled 
states extend from the plane into the pyramid of the states forming such large volumes of 
bound entanglement that include the Horodecki states as a small line. Finally, we construct 
a convex subset of the three-parameter family, a polygon, that encloses states that are 
necessarily separable. 
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